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The Impact of a Flat Plate on a Water Surface 

By J. H. G. Verhagen 1 


The phenomena occurring during the impact of a flat plate on a water surface are investi¬ 
gated both theoretically and experimentally. It is shown that, by taking the influence of 
the compressible air layer caught between the falling plate and the water surface prop¬ 
erly into account, the experimental results can be fully explained. Numerical results 
showing the effects of falling mass, drop velocity, and plate dimensions on the maximum 
impact pressure are available. 


Introduction 

During the past decades, many theoretical and ex¬ 
perimental studies have been made to clarify the basic 
properties of slamming and to correlate slamming pres¬ 
sures with the section form of a ship’s hull. 

In most theoretical treatments of the ship slamming 
problem the water is considered to be incompressible and 
nonviscous. In the mathematical idealization of the 
problem it is assumed, moreover, that the velocity of a 
specified part of the free water surface is increased in¬ 
stantaneously from zero to the downward speed of the 
falling body. The slamming pressures resulting from 
such analyses are in good agreement with experimental 
results as long as the deadrise angle of the ship’s bottom 
is not too small (see e.g. Wagner [l] 2 ). In the extreme 
case of a flat bottom, however, the theory predicts in¬ 
finitely high pressures at the instant of impact. There¬ 
fore, some buffer device should be developed to explain 
the discrepancy between theory and experiment. For 
this purpose, the compressibility of the water can be 
taken into account (von Karm&n [2], Ogilvie [3], 
Egorov [4], et al.). The initial pressure at the instant of 
impact resulting from such treatments appears to be in¬ 
dependent of the plate dimensions and equal to the 
acoustic pressure p = p/c/F, where p L is the density of 
water, c L the velocity of sound propagation in water, and 
V the downward body speed at impact. The duration of 
t he pressure peak is in the order of t = l/c h where l is 
t he half width of the flat bottom. 
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Sharov [5] considered also the influence of the elasticity 
of the construction on slamming pressures. He concluded 
from his investigation that the initial pressure is equal to 
PiCiV, independent of the elastic properties of the con¬ 
struction. 

The high pressures predicted by these theories could 
not be verified experimentally. Measurements on this 
subject show a much lower impact pressure acting, more¬ 
over, during a longer time. Fujita [6] and Tanaka [7] 
considered the influence of the air layer caught between 
the flat plate and the water surface on the slamming pres¬ 
sure. Tanaka assumed that the pressure of the viscous 
air layer between a circular plate and the water causes a 
depression in the water surface. The plate would touch 
the water annularly from its circumference to the center. 
He carried out his calculations in the same way as Wag¬ 
ner’s, modified with some of his own assumptions. 

Fujita believes that the impact of a flat-bottom body is 
due to the pressure increase caused by the rapid compres¬ 
sion of the air layer between bottom and water surface. 
He found the results obtained from such an assumption 
in good qualitative agreement with the experiments. 

In order to be able to develop a more realistic theoreti¬ 
cal model, it was deemed necessary to have a clearer pic¬ 
ture of the phenomena occurring at the moment of im¬ 
pact. Experiments with different impact velocities, 
therefore, were carried out on a rigid flat rectangular 
plate slamming on a smooth water surface. The phe¬ 
nomena were filmed with a high-speed camera (8000 
frames/sec). The pressure was recorded as a function 
of time. The pictures show that the water elevation is at 
first noticeable at the plate edges just before the plate 
touches the water surface. Next, the air layer breaks up 
into bubbles beginning at the edges and extending to the 
center of the plate with a speed in the order of magnitude 
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Fig. 1 


equal to the velocity of sound in air. With increasing 
time the bubbles drift relatively slowly outwards. 

A theory is presented to calculate the pressure as a 
function of time during impact. 

Considering the flow in the air layer under the falling 
body to be a one-dimensional compressible flow of a 
perfect gas, the flow conditions are calculated by the 
method of characteristics. The disturbance of the water 
surface due to the imposed air pressure is solved with the 
aid of potential theory. Compressibility of the water is 
neglected, as the characteristic time interval for the im¬ 
pulsive pressure appears to be large compared to l/c h 
At the moment the edges of the plate touch the upheaved 
water surface an air cushion is trapped. This volume of 
air is calculated. Assuming that the escape of air after 
that moment can be neglected and that the downward 
velocity of the mean water surface under the plate is 
independent of the place coordinate, the adiabatic com¬ 
pression of the air layer is calculated. 

Numerical results are compared with experimental re¬ 
sults. During these investigations a DTMB report has 
been published by Sheng-Lun Chuang [8] describing 
similar experiments. His results and the data reported 
here show a general agreement. 

Formulation of the Problem 

A two-dimensional plate with width 2 1 is travelling- 
downward at a velocity V(t). The flat surface of the 
plate remains parallel to the undisturbed free water sur¬ 
face, Fig. 1. The free surface is taken initially as the 
2 -axis and its shape is described by the equation 

y — Y(xjt) =0 (l) 

Basic Equations for the Air Layer 

The problem is symmetrical about the y-axis. 

The air between the plate and the water surface is 
driven away from the center to the edges of the plate. 
The pressure in the air layer increases with decreasing 
distance h(x,t) between the plate and the water surface. 
As we are especially interested in that period of time 
where the pressure increases reasonably, i.e. for values 


of h/l « 1, the flow in the air layer is considered as one¬ 
dimensional. 

Neglecting viscosity effects, the equation of motion of 
the air is 


The continuity equation becomes 

d(ph) b(puh ) _ 

hr + dx 

where 

u(Xjt) — velocity of air in horizontal direction 
p(x,t) = pressure force 
p(x,t) ~ density of air 

Assuming the compression of air to be isentropic 


Po \Po 


and using the equation 


where c a is the velocity of sound in air, the equations of 
motion and continuity can be written as: 


d x d 

dt + (u + c,,) dx 


(■■a J0/| dh\ 

h Id/ + U 5 if 


dt +(u ~ ° dx 


7 — 1 


c a jdh 
h 1 dt 


Apart from this system of partial differential equations of 
hyperbolic type, a third relation is needed for the three 
unknown variables u(x,t), c a {x,t) and h(x,t). 
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Fig. 2 


The relation between the thickness of the air layer and 
the water elevation is given by 


bh 

bt 


= V(t) + 


bY(x,t) 
bt " 


( 7 ) 


The downward velocity V(t) can be a prescribed function 
of time, or in the case of a body free falling Tinder gravity 
it can be derived from the momentum equation acting on 
the plate. With the initial condition V = 0 for t = to 
the downward velocity is 

V(t) = g(t — to) + ^ dt J i (p —* po)dx (8) 


where M is the mass of the flat-bottomed body per unit 
length and g is the acceleration of gravity. 

Since the pressure is on account of (4) and (5) related 
to the sound velocity by 



a relation between the pressure distribution in the air 
layer and the water elevation Y(x,t) suffices to determine 
the motion completely. 


important role on slamming phenomena. Therefore, the 
assumption is made that the water can be treated as an 
incompressible fluid. Assuming further nonviscous and 
irrotational flow, a velocity potential can be introduced. 
This velocity potential has to satisfy the Laplace equation 
as well as boundary conditions on the free surface. For 
t < 0 the water is undisturbed, and starting with t = 0 
the pressure over the free surface is a given function 
p (x, t). Assuming the disturbances of the water surface 
to be small, the linearized free-surface conditions can be 
used 

Y t — <p v = 0 

on the plane y — 0 

„, + !! k + £M_o (io) 

Pi 

where p t is the mass density of the liquid 

The solution of this problem is given by Wehausen and 
Laitone [9]. The complex velocity potential is 

/fe ° - c,,s 


Wave Motions Resulting from o Pressure Distribution 
Suddenly Applied to the Water Surface 

Measurements on the pressure-time history occurring 
during the impact of a flat-bottomed body on the water 
surface show that the maximum impact pressure is only 
a few kilograms per square centimeter up to impact 
velocities of several meters per second. 

The duration of the pressure impulse is long compared 
with the time required by an acoustic disturbance in the 
water to travel a distance equal to the half width of the 
flat bottom. These observations suggest that it is hardly 
likely that the compressibility of the water is playing an 


{( gk) l/t (r - t)}e~ ik(z ~®dk (11) 

where z = x + iy and y < 0 

The velocity of the water surface follows from 

(i2) 

which can be integrated to yield Y(x,t). 

The water surface elevation has to be calculated from 
t = 0 till the time t c , at which there is a first contact be¬ 
tween some point of the falling body and the water 
surface. 


DECEMBER 1967 


213 







From experimental observation it can be deduced that 
PigY(x,t) 

P(x,t) 


Due to the condition p(x,t) = 0 for \x\ > l the integration 
over £ from — oo to + °° can be reduced to the intervai 

It may be observed from equation (16) that Y(±l,t) 
will tend to infinite when bp/bx does not equal zero at the 
points £ = ±1. 


must be small compared to unity for the largest part of 
the region \x\ < l and for 0 < t < t c . On account of this 
observation and in order to simplify the calculations the 
gravity term is neglected. Equation ( 11 ) reduces then to Solution of the Aerodynamic Part of the Problem 

/m = -- 1 f dr r 

irpiJ o J - » 
and equation ( 12 ) becomes 

rV Y(x, l) = — lim 1 Im 

uo V TO TTPl 

With the condition that p( 
tegration by parts yields 

bY(x,l) __ 1 

bt “ 


TTPl 


W(Z,t) ,, , 

~~zzj * or y < o 

(13) 

We have to solve the following system of partial dif¬ 
ferential equations for the unknown functions u(x,t ), 
c a (x,t), and h(x,t). 


(14) 

(d/ + (” ■*" '") dr} (“ + 7 1 ) 


x,l) is zero at infinity, 

an in- 

c a j dh 

~~ ~ h \dl + 

bh) 

%,} < 1? ) 

bp 

fur r * 

[Jo J -co x — £ 

(15) 

{1 + (M “ rJ dtj (“ ~ 7 i) 



Thus, the elevation of the water surface becomes 


c a jbh bh 
h 1 d£ U bx 


(18) 


bp 


Y(x,t) 


1 r r r a 

J7 I dr I dr' I -— d£ (16) 

T Pl Jo Jo J -oo X — £ K 7 


A third relation between h(x y t) and c a (x,t) can be ob¬ 
tained by integrating relation ( 7 ). With the use of equa¬ 
tions (16) and ( 9 ) the result is 
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with p(dtzl } t) = p 0 , is 



(19) 


where ho is the position of the plate at time t — 0, while 
t ~ 0 is a rather arbitrary starting point of the calculation. 

Initial and Boundary Conditions 

To determine the solution of the foregoing hyperbolic 
differential equations, prescribed values of u, c a , and h 
along a curve in the x-t plane are needed. 

Due to the formulation of the problem it must be as¬ 
sumed that the validity of the derived equations is re¬ 
stricted to that region of the x-t plane for which \x\ < L 
The time interval is also restricted by the condition that 
h/l must be small compared to unity. 

At first glance, it seems obvious to make the following 
simple choice of the initial conditions on the line segment 
t = 0, \x\ < l: 

h(x,o) — ho , u(x,o) = 0, c a (x,o) = r 0 


where c 0 is the acoustic speed of sound in air. The last 
condition implies that the pressure under the plate is 
equal to the undisturbed pressure in the surrounding 
medium. However, there are some objections against 
these particular conditions. First, the question can be 
raised whether the final result will be affected by the arbi¬ 
trary choice of the initial conditions. Secondly, these 
conditions lead to the appearance of a shock wave in the 
air layer, travelling from the edges of the plate to its 
center, with the acoustic speed Co. The pressure dis¬ 
tribution on the water surface contains now a disconti¬ 
nuity due to the pressure jump across the shock wave, 
which results in an infinite water surface elevation accord¬ 
ing to equation (16). 

To avoid mathematical difficulties and to get a more 
realistic picture of the physical phenomena, another 
choice of the initial conditions is made. 

If the velocity of the air under the plate is small com¬ 
pared with the velocity of sound Co, it may be assumed 
that the air is incompressible and the water surface un¬ 
disturbed. The air velocity, which follows from the con¬ 
tinuity equation 


dh , l du A 

77 + h - . = 0 

dt dx 


with u(o,t) = 0, is 


u(x,t) 


x dh 
h dt 


( 20 ) 


and the pressure distribution, from 


po dx 


du 

dt 


-T u 


du 

dx 


l 2 _ x 2 

p(x,t) - Po = po — 



where 


hd/h) 
2 dt 2 } 


h(t) = h 0 



( 21 ) 


( 22 ) 


With the downward velocity V(t) of the plate given, the 
flow conditions can be calculated. 

The choice of ho is restricted by the conditions that 
ho/1 « 1 and u(l, t)/c 0 <K 1, so the initial conditions for 
— l< x < + 1 , t = 0 can be written as 


— «T « 1; Y(x,o) = 0 (23) 

Co l 



With these initial conditions no shock wave in the air 
layer will appear. The initial excess pressure under the 
plate is small enough to justify the assumption that the 
water surface is undisturbed. This implies that the 
effect of the choice of h 0 on the final result can be neg¬ 
lected as well. ( N.B . h 0 is the position of the plate at 
the starting point of the numerical calculation; it has 
nothing to do with the total drop height of the plate.) 

The solution of equations (17-19) can be determined 
now in the region I of the x-t plane bounded by two 
characteristics, defined by dx/dt = u ± c a and the line 
segment t = 0, — l < x < +1 (see Fig. 2). 

In order to extend the solution to larger values of time, 
boundary conditions for x = ± l are needed. 

It seems obvious to assume that p(±l, t) = po and 
hence c(±Z, t) = c 0 as long as the escape velocity of the 
air at the edges of the plate does not exceed the acoustic 
sound velocity c 0 . 

Let t' be the time at which the escape velocity of the 
air attains the value Cq. For t > t' the escape velocity at 
the plate edges will be equal to the local sound velocity 
c a (l,t ). So the boundary conditions are: 


for 0 < t < t' 

CM ) = C a (-l,t) = Co (26) 
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Fig. 4 


for t > V 

u(l,l) = c a (l,t) 

-u(-l,t ) = c.(-M) (27) 

where the value of t' is defined by 

u(l,l') = Co (28) 

To avoid inconvenient singularities in the mathemati¬ 
cal description of the problem, the conditions at the 
edges of the plate need a slight modification, as will be 
discussed in the next section. 

Flow Condition at the Edges of the Plate 

In the previous sections use was made of the statement 
that the pressure on the linearized free water surface 
immediately outside the plate edges is equal to the un¬ 
disturbed pressure prevailing in the surrounding medium, 
so 

p(x,y,t) v -o = po for \x\ > l 

When the undisturbed water surface is replaced by a 
rigid bottom and the pressure under the plate is calculated 
with the use of the boundary conditions (26-28), the re¬ 
sult will be as given in Fig. 3. The first derivative of the 
pressure distribution in the ^-direction will in the first 
case show a jump and in the second case a delta function 
at the edges of the plate. According to equation (16) the 
water elevation should have a logarithmic singularity lor 
t < V and a l/x singularity for t > l '. Actually the pres¬ 
sure distribution will be continuous at the edges of the 
plate. The used approximation of one-dimensional flow 
in the air layer is too rough at the end points. T he only 
conclusion that can be drawn from the foregoing is that 
the acceleration of the water surface at the plate edges is 
much more important for l> V than for smaller values of 
time. Consider the air flow of the jet issuing from the 
slit between the edges of the plate and the water surface 
more closely. The Reynolds number uh/v with which the 
plane jet enters into the surrounding medium is in prac¬ 
tical cases sufficiently large to ensure the jet to be fully 
turbulent. Suppose the water surface remains undis¬ 


turbed. The structure of the jet may be sketched as in 
Fig. 4. Near the slit, the jet boundary is diffused into a 
highly turbulent plane “mixing zone” with wedge-shaped 
profile. This gradually “eats” into the undisturbed 
central core of the jet. The angle of spread of the mixing 
zone ranges from 25-30 deg, and causes the nonturbulent 
cone of the jet to be eaten up after four or five times the 
slit opening [10]. At the boundaries of the jet the pres¬ 
sure equals p 0 . Inside the cone on the line y — 0 the 
pressure will gradually decrease in the x-direction to the 
value p 0 . 

In view of this explanation it is assumed that the pres¬ 
sure curve as a function of x as shown in Fig. 3 must be 
smoothed at the plate edges over a region proportional 
to h, say ah where a is of the order of unity. 

Owing to this revised pressure distribution, the singu¬ 
larities in the water elevation at the plate edges are trans¬ 
formed into finite peak values. From equation (16) to¬ 
gether with the conditions as stated in the foregoing, the 
maximum surface elevation can be approximated by 


7(±M) = — f dr Pin 

7Tp/ JO 


ah\ dp(l,o_y) 


r(*v> - W> + „*,£* JV 


for t > t! 


With increasing time, the elevation Y(l,t) increases as 
well. A moment will be reached at which Y(l,t) will be 


equal to h Q - J Vdt. At that moment the air layer 

under the plate is cut off from the surrounding 
medium; and the air flow issuing from the slit will be re¬ 
duced suddenly to zero. 

It is highly probable that a shock wave is formed in the 
air layer as soon as the free escape of the air is cut off by 
the disturbance of the water surface at the plate edges. 
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5 Formation of bubble layer at moment of impact. The views are taken at 8000 frames /sec from beneath water surface 


These shock waves travel from the edges to the center of 
the plate with the local sound velocity. The pressure 
over the shock wave is discontinuous, resulting in such a 
large water surface disturbance (equation 16) that the 
air layer is broken up into bubbles after the passage of 
the shock front. This is supported by experimental 
evidence. 

Fig. 5 shows a series of pictures taken with a high¬ 
speed camera at 8000 frames/sec. The lower side of the 
plate is filmed from beneath the water surface. The 
velocity at which the bubble formation is expanding as 
taken from these pictures is about 350 m/sec. The cal¬ 
culated value of the mean sound velocity for this case is 
Cmean = 1-08 Co & 360 m/sec. The plate touches the 
water surface as soon as 

Y(l,t) = h 0 - J'vdt. 

The description of the slamming problem after that time 
needs another mathematical idealization, which will be 
described later in the paper. 


Numerical Solution of Hyperbolic Partial 
Differential Equations 

The quasi-linear system of equations for the dependent 
variables u and c a as functions of the independent vari¬ 
ables x and t as stated in (17) and (18) has to be solved. 

With the use of equation (19) the right-hand sides of 
equations (17) and (18) can be expressed in terms of u, c a , 
x, and t. This expression appears independent of the 
derivatives of u and c a . So there are two characteristic 
directions 


dx 

dt 


u + c a 


dx 

dt 


u — c a 


(30) 


The characteristic equations are 

2 c a fbh dh\ 1 

du + Tzi dc * + h ( dt +u Sx) dt = 0 

(31) 

2 c a fdh dh\ 

du ~7=i dc °-h\bt +u zx) dt=:0 

where h is given by equation (19). 

Because the boundary and initial conditions of the 


problem are given on lines parallel to the axes of the x-L 
plane, the numerical method of specified time intervals is 
preferred instead of using a grid of characteristics. 

The method of specified time interval starts with the 
assumption that u and c a are known functions of x at a 
specified time t = t , either as given initial conditions or as 
the result of a previous stage of the calculation. Using 
finite-difference approximations to the two partial dif¬ 
ferential equations along the characteristic directions, the 
values of u and c a can be obtained at the time t = t + At. 
The computational procedure is described by Mary Lister 
[11] and will not be repeated here. The boundary con¬ 
ditions on the lines x = ±Z are stated in (26)-(28) and 
the initial conditions for —Z < x < l in (23)-(25). The 
correction in the pressure distribution at the edges of the 
plate is described in the previous section and is applied 
only to the calculation of the surface elevation and hence 
to the local distance h(x,t) between plate and water. 

The factor a, which is a measure for the region over 
which the pressure curve is smoothed at the end points, 
appears in these expressions as a quantity of the order 
unity whose exact numerical value is unknown. In the 
numerical computation the value of the factor a was 
varied. Its influence on the final result turned out to be 
rather small. The reason for this is connected with the 
following important conclusion which could be drawn 
from the numerical computation. The moment at which 
the air layer under the plate is locked appears to coincide 
almost with the moment that the escape velocity of the 
air attains the acoustic sound velocity. The influence 
of the numerical value of a on this phenomenon is rather 
small, which can be explained as follows. When the 
escape velocity of the air is subsonic the acceleration of 
the water surface under the plate edges is proportional 
to In ( ah/l ). After the acoustic sound velocity is at¬ 
tained the proportionality changes into ( ah/l)~' 1 [see 
equation (29) ]. Because ah/l «< 1 at the moment that 
u(ld) = c 0 , the elevation of the water surface is much 
more important in the latter case. 


The elevation Y(Id) attains the value ho ~ 


/: 


Vdt 


just after the critical escape velocity is attained. A 
variation in the choice of a causes only a slight shift in 
this result. So the fact that the exact value of a is un¬ 
known is not so important for the practical solution of the 
problem. 
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Fig. 8 Flat-plate model 


Installation of test model 


Compression of the Trapped Air Cushion 


At last the description is given of a mathematical model in the time intervals under consideration th 
for the air cushion, which elucidates its action as a the plate will be small in general compared to 
buffer device. When the two-dimensional plate catches load, 
the water surface at its edges, the mean thickness of the ^y 

trapped air layer is hi, and the mean pressure is equal to M ^ = —2 l(p — p 0 ) 

pi. It is assumed that from this moment the downward 

velocity of the water surface under the plate is only a The same equation for the virtual mass yields 

function of time and independent of the ^-coordinate. ^ w 

Further, the flow of air in the ^-direction is neglected. m ^ = 2 l(p — p 0 ) 

The pressure in the air or bubble layer is thus place- 

independent. Due to the assumption that the velocity The air compression will be isentropic, so 
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Fig. 9 Impact pressure record for test with drop height of 4 and 40 cm respectively 


Using these equations, the following relation for the 
air layer thickness is obtained: 

d 2 h div dV m + M 

dt 2 dt dt mM ^ 

The initial conditions are 

t = h] V = Vi) w = wi (37) 

The differential equation (36) is nonlinear. 

By multiplying both sides of the equation by 2 dh/dt 
and carrying out an integration over t, a linear differen¬ 
tial equation is obtained. 


h 


Po 

PiJ 


(36) 



m + M 
mM 


4/pi 



Vo 

ViJ 


dh 7 

n di 


+ (wi - Fx ) 2 

Evaluating the right-hand side gives 


Using equations (32) and (35) and after rearranging, this 
transforms into 

d(y/pi) 
d(c 0 t/l ) 



+ M 


M 


Sylpg I" Pi 
7T hpi \_po(y — 1) 



+ 



(39) 


1 / dh\ 2 m + M 4/Aipo Pi 

c 0 2 \dt) M mco 2 _ p 0 (y — 1 ) 



With the initial condition p = pi for t = ti the relation 
between pressure and time can be solved numerically. 
/on\ The impulsive pressure attains a maximum when dp/dt 
^ =0. The magnitude of p = p m&K follows from: 
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it pi yl M 



(40) 


The initial conditions of the trapped air layer hi/l, pi/po, 
Wx/cq, and Vi/co can be obtained from the calculation 
described in the previous sections. 

The positive sign in equation (39) corresponds to a 
pressure increase and has to be taken before the maxi¬ 
mum pressure is reached. Afterwards the negative sign 
must be used. Moreover, the right-hand side of (39) is 
only a function of p, for given initial conditions. There¬ 
fore, the pressure-time curve is symmetric around the 
time at which the maximum pressure is attained. This 
conclusion is in agreement with the measured pressure¬ 
time history curve. 


Description of Experimental Setup 

A cast aluminium plate 40 by 40 by 2 cm with bar 
stiffeners as shown in Fig. 7 is dropped from various 
heights up to 50 cm on a smooth water surface. The 
total weight for the test was 8 kg. The tests were carried 
out in an aquarium 100 by 40.5 by 40 cm with glass walls 
and bottom. 

The support consisted of a steel frame which guided a 
vertical shaft, at the end of which the flat plate was 
attached. To simulate two-dimensional flow conditions 
the plate was provided with perspex vertical strips at two 
parallel edges; see Fig. 8. A pressure transducer of 
barium-titanaat crystal was fitted at the center of the 
plate. 

By means of a photoelectric cell the point at which the 
plate had fallen to a few mm above the water could be 
marked. At that moment the time base of the oscillo¬ 
scope was started and the pressure-time history photo¬ 
graphed. Fig. 9 shows typical pressure-time history 
curves during a 4-cm and a 40-cm free fall. The maxi¬ 
mum values of the impact pressure are plotted as a func¬ 
tion of the drop height in Fig. 10. 


Numerical Results and Experimental Data 

The effect of the following variables on the impact 
pressure is numerically evaluated, where: 

h — drop height 
l = half width of the flat bottom 
M = mass of the free-falling body 

Normal atmospheric conditions are chosen, namely 
g = 9.81 m/sec 2 

pi = 10 3 kg/m 3 

p a = 1.29 kg/m 3 ] 

pa = 10 5 N/m 2 / —► c a = 330 m/sec 
y = 1.40 ) 

Fig. 10 presents the results for two bodies of different 
masses M/pil 2 = 0.5 and 1, respectively. The maximum 


impact pressure is plotted as a function of the drop 
height. The half width of the flat bottom is Z = 0.2 m. 
It is seen that p max is approximately directly proportional 
to the drop height or to the square of the downward body 
speed. This proportionality appears to be valid for small 
values of M/pil 2 . With increasing mass, however, there 
is a more or less linear relationship between p max and F, 
which can be concluded from the following considerations. 

The most important pressure rise under the body 
occurs after the moment the escape velocity of the air 
reaches the velocity of sound. As follows from mass con¬ 
tinuity, the thickness of the air layer at that moment is 
of the order h/l = V/c a . The body travels down over a 
distance h during a time of the order A t = h/V = l/c a . 
If M/pil 2 )$> 1, the body velocity is hardly decelerated by 
the influence of the air cushion, so the added mass of the 
body has to accelerate to the velocity V during the same 
time interval At. 

The total impulse produced in the system will be ap¬ 
proximately equal to the change in momentum of the 
added mass m alone 

r*t o-j- At 

I = 21 l pdt mV 

J to 

Hence the maximum impact pressure p max is proportional 
to 


Pmax 


mV 

~lAt 


PiVc 


a 


which is linear with the downward body speed. 

Fig. 11 indicates that the maximum pressure is inde¬ 
pendent of the mass of the falling body when M/pil 2 » 1. 
The diagram is calculated for a constant drop height h = 
0.2 m (Z = 0.2 m). This result is qualitatively in 
agreement with experimental data on the impact of cir¬ 
cular plates of varying mass on the water surface, ob¬ 
tained from Y. Fujita [6]. 

Fig. 12 shows the pressure during impact as a function 
of time for two different drop heights. The numerical 
results are compared with experimental data in Figs. 10 
and 12 and the agreement is excellent 

Finally, Fig. 13 shows the influence of plate dimension 
on the maximum impact pressure. The value of M/pil 2 
= 1 and the drop height h = 0.2 m. It appears that the 
maximum impact pressure is nearly independent of the 
plate width for equal values of the ratio M/pil 2 . 

Hence, to obtain the same impact pressure of two 
bodies falling from an equal distance above the water 
surface, the mass of the body per unit bottom surface 
must be directly proportional to the two-dimensional 
body width. The larger body must be relatively heavier 
than the smaller one. 


Discussion of Validity 

The most important assumption lies in neglecting the 
compressibility of the water and the elasticity of the 
body. Compressibility effects are neglected because the 
events we are interested in are expected to happen in 
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times of the order required by an acoustic wave in the air 
to travel over a distance l , i.e., At = l/c a , which is large 
compared with l/c h The experiments carried out here 
indicate that this assumption is fully justified. 

The excellent agreement between theory and experi¬ 


ment shown in this paper may suggest that the subject of 
flat-bottomed-body slamming has now been exhaustively 
described. The experiments, however, are limited to 
small values of the mass of the body compared with the 
added mass. The author feels that the theory enunciated 
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half width of Flat bottom 

Fig. 13 Effect of plate width on maximum pressure 


here will be valid only for not too high values of M/pf 2 . 

With increasing mass of the body the maximum impact 
pressure increases. When M/pf 2 1 the pressure-time 
curve starts approximately like (c a /l)(t — r) 1A , where r 
is the time at which the plate should reach the water 
surface in the absence of the air layer. With increasing 
value of the maximum pressure the gradient of the pres¬ 
sure-time curve becomes steeper. An important pressure 
variation can occur in a time interval small compared 
with c a t/l, such that the time-dependent term in the 
wave equation can no longer be neglected. That is, the 
compressibility effects of the water cannot be neglected 
when the maximum pressure increases. Also, the elas¬ 
ticity of the construction in that case must be taken into 
account. 

A comparison of experimental data on a 20 in. by 26.5 
in. flat-bottom model M/piL 2 & 2, obtained fromSheng- 
Lun Chuang [8], shows a considerably lower maximum 
impact pressure than predicted by theory; see Fig. 10. 
The discrepancy may be due to the neglected effect of 
water compressibility mentioned before, as well as to the 
possibility that the experimental model (a stiffened steel 
box) may not be sufficiently rigid to withstand the im¬ 
pulsive pressure without deflection. An indication of this 
J possibility could be that his measured maximum impact 
pressure appears somewhat higher at the edges of the 
model than at the center, which is not to be expected 
[6]. Further, the impact acceleration history curve of the 
model shows a discrepancy in shape with the impact 
pressure history, which is in contradiction to Newton’s 
second law for a rigid body. The time duration of the 
impact pressure is in agreement with the author’s ex¬ 


periments, which should justify also in this case our 
neglect of water-compressibility effects. 
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